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OBSTACLES  AS  PROBES  OF  THE  BLAST  WAVE 
INTERIOR  IN  THE  NRL  LASER/HANE 
SIMULATION  EXPERIMENT 


I.  Introduction 


Over  the  past  several  years  the  Plasma  Physics  Division  at  the  Naval 
Research  Laboratory  has  developed  an  experimental  apparatus  which  can  simulate 
certain  physical  processes  characteristic  of  a  high  altitude  nuclear  event 
(HANE).1  This  program  is  under  the  auspices  of  the  Division  of  Atmospheric 
Effects  of  the  Defense  Nuclear  Agency.  It  has  the  primary  goal  of  a  physical 
understanding  and  reliable  prediction  of  the  degradation  of  radar  and 
communication  systems  through  a  HANE-disturbed  atmosphere. 

The  first  series  of  experiments  (prior  to  January  1985)  were  initiated  by 
irradiating  a  thin  target  with  a  single  neodymiun  laser  beam.  Dark-field 
shadowgraphy  for  the  case  of  a  high  pressure  (2  1  Torr))  ambient  gas  showed 
that  the  rapid  heating  and  evaporation  of  the  target  material  was  followed  by 
an  expanding  shock  wave.2,3.  Thus  in  a  real  HANE,  as  well  as  in  the 
experiment,  the  disturbed  region  was  bounded  on  the  exterior  by  a  strong  shock 
front.  However,  the  velocity  and  temperature  within  the  disturbed  region  of 
the  experimentally  produced  explosion  has  not  been  discerned.  We  will  refer 
to  this  interior  region  as  the  cavity  because  the  electron  density  has  been 
found  to  be  much  lower  in  the  interior  than  at  the  shock  front.1*  Clearly  the 
degree  to  which  the  experiment  can  simulate  a  HANE  and  the  applicability  of 
the  laboratory  results  for  predicting  nuclear  environments  depends  on  a 
knowledge  of  the  flow  within  the  cavity.  The  present  report  is  addressed  to 
this  problem  of  the  cavity  dynamics  and  will  propose  a  new  experimental 
diagnostic  to  directly  measure  the  Mach  number  in  the  cavity. 

As  background  let  us  now  discuss  the  present  evidence  for  the  flow 
conditions  within  this  cavity.  From  the  analysis  of  the  experimental  data  by 
Ripin  et  al.^,  it  appears  that  the  shock  front  and  attendent  high  density 
shell  expand  like 
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where  Rg  is  the  shock  radius  in  cm,  E  is  the  laser  energy  in  Joules,  pQ  is  the 
ambient  pressure  in  Torrs,  u  is  the  mean  molecular  weight,  and  t  is  the  time 
in  nsec.  This  formula  has  exactly  the  same  functional  dependencies  as  the 
spherical,  adiabatic,  expansion  of  a  Taylor-Sedov  self-similar  blast  wave.® 
In  this  ideal  theoretical  model  and  for  a  ratio  of  specif c  heat  (Y)  between  1 
and  3,  the  density  decreases  and  the  temperature  increases  toward  the  center 
in  such  a  way  that  the  pressure  is  roughly  constant,  while  the  velocity 
decreases  nearly  linearly  with  radius  inward.  If  electron  thermal  conduction 
is  important  the  adiabatic  assumption  can  be  replaced  by  an  isothermal  one  in 
the  extreme  limit.  In  this  case  the  expansion  is  still  self-similar,  with 
again  the  functional  form  of  eqn  (1),  but  the  interior  temperature  is  that  of 
the  shock  front  and  the  density  and  velocity  decrease  toward  the  center.® 

The  condition  for  applicability  of  the  self-similar,  Taylor-Sedov 
solution  to  an  explosion  has  generally  been  that  the  initial  conditions  are 
forgotten.  This  can  be  translated  into  the  statement  that  the  mass  3wept  up 
by  the  shock  front  (M  -  4irp  R  ®/3,  where  p  is  the  ambient  density)  be  much 

SW  0  3  0 

larger  than  the  debris  mass  (Md).  Detailed  numerical  simulations  by  Gull' 
show  that  by  the  time  Mgw/Md  -  10  the  outer  region  of  the  cavity  resembles  the 

Taylor-Sedov  blast  wave  solution  while  the  inner  region  does  so  by  the  time 

M3w/Md  ~  50*  111  the  experiment  a  typical  target  mass  is  4ugr  but 

only  -  0.3ugr  are  ablated  to  become  debris.  Using  the  relation 

■3  _c  u(Np) 

Po(gr/cnr)  «  1.6  x  10  pQ(T)  — “ —  (2) 

at  standard  temperature,  one  finds  for  pQ  -  5  Torr  that  M3W/Md  *  10 

at  R  -  0.45  cm,  and  MaiJ/MH  -  50  at  R  -  0.76  cm.  For  reference  at  23  Joules 

of  laser  energy  the  corresponding  times  are  12  nsec  and  44  nsec, 

respectively.  Shadowgraph  observations  typically  lie  within  55  to  155  nsec. 

Given  the  above  two  facts  for  the  high  pressure  runs;  vis.,  (i)  that  the 

blast  wave  radius  follows  the  Taylor-Sedov  relation,  and  (ii)  that  the 

ratio  M  /M .  >>  1  at  R  -  1  cm,  one  is  led  to  conclude  that  the  flow  within 
sw  d  s 

the  cavity  follows  the  adiabatic  (or  isothermal)  Taylor-Sedov  solution. 
Consequently,  the  cavity  should  be  at  least  as  hot  as  the  shock  front. 
However,  recent  numerical  calculations  of  the  experiment  by  Stellingwerf0 
predict  a  quite  different  picture.  Using  the  same  parameters  mentioned 

above  (pQ  -  5  Torr  of  N 2  and  EQ  -23  Joules),  he  found  that  the  temperature 


decreases  from  the  swept-up  shell  back  into  the  interior,  i.e.,  a  cold  cavity, 
even  as  late  as  64  nsec.  He  finds  no  evidence  that  the  solution  is  evolving 
toward  the  Taylor-Sedov  blast  wave,  although  the  blast  wave  radius  does  agree 
with  the  t dependence  of  eqn.  (1).  The  reason  for  the  discrepancy  from  the 
Taylor-Sedov  blast  wave  is  not  obvious  from  the  calculation.  It  is  noteworthy 
that  the  reflected  shock,  which  forms  at  one  equal  mass  radius,  fails  to  move 
back  into  the  cavity  and  heat  it  as  would  be  expected  in  the  approach  to  the 
Taylor-Sedov  solution.  In  any  case  it  is  clear  that  the  prediction  of  a  cold 
cavity  is  in  complete  variance  with  the  Taylor-Sedov  model. 

Unfortunately,  there  are  no  direct  experimental  observations  on  the 
temperature  of  the  cavity  to  distinguish  between  the  two  above  models. 
Spectroscopy  could  in  principle  measure  the  thermodynamic  state  in  the  cavity 
but  the  measurements  would  be  severely  hampered  by  the  strong  continuum  from 
the  swept-up  shell.  In  one  shot  (#14197)  a  large  wedge  was  placed  so  that  it 
protruded  into  the  cavity  by  155  nsec.  The  original  intent  was  to  see  if  the 
aneurism  was  affected  (it  wasn’t).  The  resulting  non-statlonary  oblique 
shocks  seen  in  the  dark-field  shadowgraph  are  interesting,  but  the 
configuration  of  the  shocks  and  wedge  was  so  complex  that  any  conclusion  based 
on  the  shock  geometry  is  highly  tentative. 

Taking  off  from  this  idea,  though,  we  suggest  the  use  of  small  obstacles, 

such  as  spheres  or  conical  bullets,  to  directly  probe  the  cavity.  Consider 

the  sequential  Interaction  of  a  moving  shock  front  with  a  sphere  fixed  in  the 

lab  frame  as  depicted  in  Fig.  1 .  As  the  main  shock  wave  engulfs  the  obstacle 

q 

a  reflected  shock  is  formed  and  meets  the  main  shock  at  the  triple  points.7 
Further  into  the  sequence  the  reflected  shock  expands  away  from  the  obstacle 
and  the  main  shock  reforms  downstream.  At  this  point  the  obstacle  is  subject 
to  the  postshock  flow,  which  relative  to  the  main  shock  is  subsonic,  but  may 
be  subsonic  or  supersonic  relative  to  the  obstacle.  Let  M.  be  the  Mach  number 
of  the  postshock  flow  relative  to  the  obstacle.  If  M1  is  <  1  the  reflected 
shock  continues  to  move  upstream  of  the  obstacle  and  degrades  into  a  sound 
wave.  If  M1  >  1,  a  standoff  bow  shock  is  formed,  as  shown  in  the  last 
sequence  of  Fig.  1.  In  this  case,  classical  gas  dynamics  provides  a 
quantitative  relation  between  the  standoff  distance  and  the  Mach  number  M1 . 
Thus  once  an  obstacle  is  engulfed  by  the  expanding  cavity,  the  absence  or 
presence  of  a  bow  shock  and  its  standoff  distance  can  indicate  the  local  Mach 
number.  To  go  from  here  to  a  local  temperature,  an  assumption  of  the  velocity 


and  ratio  of  specific  heat  is  required.  Since  both  the  Taylor-Sedov  solution 
and  Stellingwerf's  calculation  have  a  nearly  linearly  velocity  profile,  the 
use  of  a  linear  velocity  law  for  the  experimental  data  is  reasonable. 
Actually,  the  procedure  could  be  reversed  to  provide  a  further  test  on  the 
theoretical  and  numerical  models;  from  the  calculated  Mach  number  of  a 
numerical  model  the  standoff  distance  of  a  bow  shock  from  a  small  obstacle 
could  be  predicted  and  compared  with  experimental  observations. 

We  present  a  cursory  review  of  the  gas  dynamic  theory  for  standoff  bow 
shocks  in  section  II.  The  main  results  are  presented  in  Figures  5a  and  5b. 
The  applicability  of  the  theory  to  the  experiment  is  discussed  in  section 
III.  There  we  first  list  the  assumptions  behind  the  simple  theory  and  the 
complications  in  the  experiment  Then  we  note  that  it  would  be  useful  to  study 
the  temporal  evolution  of  the  standoff  distance  as  a  gross  analysis  of  the 
cavity  dynamics.  This  approach  is  exemplified  by  discussing  the  problem  for 
the  adiabatic  and  lsothemal  self-similar  blast  waves,  for  Stellingwerf ' s 
model,  and  for  a  recent  numerical  model  developed  at  NRL.  Section  IV  contains 
a  summary  with  the  major  conclusions. 

II .  ANALYSIS 

The  theoretical  analysis  to  determine  the  standoff  distance  of  a  bow 
shock  cannot  be  reduced  to  a  single  equation,  but  is  instead  composed  of 
several  separate  steps.  We  will  not  present  each  step  in  complete  detail  for 
the  basic  procedure  is  outlined  in  several  sections  of  A.H.  Shapiro's  text  on 
gas  dynamics.10  Our  approach  will  be  to  gather  the  separate  steps  into  an 
organized  sequence,  derive  or  quote  the  main  equations,  and  discuss  our  method 
of  solution.  Shapiro's  results  are  obtained  for  a  ratio  of  specific 
heat  Y  equal  to  1  .H  since  the  problem  was  first  solved  in  aeronautics.  We 
will  present  results  for  Y  -  1.2  and  5/3  to  explicitly  demonstrate  the 
dependency  on  Y.  This  range  is  considered  for  the  Y  of  the  cavity  is  unknown: 
at  the  main  shock  where  there  are  rapid  chemical  reactions  Y  -  1 .2,  while  if 
there  are  no  Internal  degrees  of  freedom  Y  -  5/3. 

A.  Maximum  turning  angle  at  an  oblique  shock 

The  first  step  is  to  determine  the  maximum  turning  angle  a  streamline  of 
incident  Mach  number  M^  can  undergo  upon  passing  through  an  oblique  shock. 


The  jump  relations  for  a  perfect  gas  across  an  oblique  shock  can  most  easily 
be  obtained  from  the  normal  shock  relations  by  letting  V1  -*■  V1  sin  0  and 
V2  -*•  V2  sin  ($  -  x).  »  where  the  geometry  and  nomenclature  are  contained  in 
Figure  2.  This  amounts  to  rewriting  the  standard  relations  in  terms  of  the 
normal  velocity  components.  One  then  finds 


V.  sin* 


(Y  ♦  1)M,2sin2d> 


>1  V^sinC  )  (Y-DM^sin2#  +  2 


2YM12sin24i 
Y  +  1 


Y  -  1 

Y  +■  1  ’ 


?  ?  2  ♦  (Y  -  1)  M  2sin2<fr 

Mp^sin^Co  -  X)  -  - p - 5 - 1 - : 

2YM  sin  *  -  (Y  -  1) 


and  the  continuity  of  the  tangential  velocity  component  is 


V1  CO30  -  7 ^  cos  (<fr  -  x)* 


By  combining  relations  (3)  and  (6)  one  readily  finds 


cot  (<fr  -  x)  »  cot 


(Y  ♦  1)  M^sin2* 

(Y  -  1)  M^sin2*  ♦  2 


from  which 


(Y  ♦  1)  M. 

cot  x  *  tan  0  [ - = - = - l]  (8) 

2(M1  sin  ♦  -  1) 

A  relation  which  will  be  needed  later  is  obtained  from  eqns.  (5)  and  (7) 


2  ♦  (Y  -  1)  M, ‘ 


2  2 
2M,  cos  ♦ 


2YM1  sin  ♦  -  (Y  -  1)  2  +  (Y  -  1)  M^  sin  0 

From  eqn  (8)  one  finds,  in  a  plot  of  x  versus  0,  that  a  line  of  constant 
M1  begins  at  (x,  ♦)  -  CO,  arcsin  ( 1  /M^ )  ]  corresponding  to  a  Mach  wave, 
increases  to  a  maximum  x_-  at  0  ,  and  then  ends  at  (x,  0)  -  (0,  ir/2) 

uiaX  ui3X 

corresponding  to  normal  shock.  The  value  of  x  can  be  found  by 
differentiating  eqn  (8)  with  respect  to  0  and  solving  for  0  : 


[(Y *  -~M1 2  ~  4)2  *  4y(l  *  — — 2~M  f)]1/2}]1/2  (10) 


The  value  of  x__v  is  then  given  by  eqn.  (8).  For  a  given  free  stream  Mach 

HlaX 

number  M« ,  1/M1  <  sin  $  S  sin  $  is  the  range  of  allowable  *  for  the  so- 
i  i  max 

called  weak  shock  solution.  The  strong  shock  solution,  wherein  $  >  *  ,  are 

T  max 

not  observed  in  the  laboratory  but  are  indicative  of  shock  detachment  from  a 
cone  or  wedge. 

B.  Supersonic  flow  past  a  cone. 


We  next  need  to  solve  for  the  flow  past  a  cone  in  a  normal  incidence 

supersonic  stream  of  Mach  number  M1 .  Given  -that  the  flow  between  the 

standing,  attached,  conical  shock  and  the  cone's  surface  is  steady,  adiabatic, 
and  irrotatlonal ,  the  gas  dynamic  equations  can  be  reduced  to  a  pair  of 

ordinary  differential  equations.  The  flow  is  conical,  i.e.,  dependent  only  on 
the  angle  9  of  Figure  2.  The  derivation  is  given  by  Shapiro10  on  page  654  in 
a  spherical  coordinate  system  with  velocity  components  Vp  and  V0  at  an 
arbitrary  point  P.  In  our  notation  the  equations  become 


d  U  Y  +  1  _  Y  -  1 

- [ -  (4°12 - 

.a2  L  2  ld0J  2 
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,.2,  du 
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vr/vmax  and  vmax  i3  the  maxlmum  velocity  for  a  given  stagnation 


temperature.  It  is  related  to  the  adiabatic  sound  speed  through 


c2  -  l-L-ifV2  -  v2) 
G  2  ' max  v 
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and  hence,  in  terms  of  a  local  Mach  number 


7  (Y  -  1 )M 

2  2' 

7  2  ♦  (Y  -  1)M 

max 

Equation  (11)  was  first  derived  by  Taylor  and  Maccoll11  in  1933  and  they 

found  excellent  agreement  between  their  solution  and  laboratory  experiments. 

Their  method  of  solution,  described  in  Shapiro,  begins  at  the  cone  surface  and 

is  cumbersome  given  modern  day  computers.  Instead  we  employ  the  following 

method.  For  a  chosen  Y,  M1  ,  and  $,  eqns,  (8)  and  (9)  give  the  turning 

angle  x  and  the  immediate  postshock  Mach  number  M2.  From  eqn.  (13)  one  gets 
2  2 

V2  ^Vmax  and  ttle  velocity  components  at  the  shock  front,  where  9  >»  4>,  are 

Vr  lV2l 

T. -  -  t. -  COS($  ~  X)  -  U, 


7 -  -  -  y -  Sin  (4i  -  X)  • 

max  max 

Finally  from  eqn  (12)  one  has  dU/d  9,  and  eqn  (11)  can  be  integrated  starting 
at  the  shock  in  the  decreasing  9-direction  till  dU/d9  «  0.  Since  this  last 
relation  means  7-0,  it  represents  the  inner  boundary  condition  on  the  cone 
surface  and  the  corresponding  angle  9  is  the  half  angle  5  of  the  cone. 

The  results  for  Y  -  1 .2  and  5/3  are  presented  in  Figures  3a  and  3b, 
respectively.  The  upper  dashed  line  in  each  graph  is  the  maximum  shock 
angle  4  for  a  given  free  stream  Mach  number  M. .  Along  this  line  the 

ludlX  ' 

starting  conditions  for  the  integration  of  eqn.  (11)  correspond  to  the 


solution  for  the  maximum  turning  angle  X 


discussed  in  the  previous 


subsection.  Consider  a  cone  at  a  fixed  half  angle  5,  i.e.,  a  solid  line.  As 
M1  decreases  the  shock  angle  <p  increases  until  is  reached.  Further 
reduction  of  M1  leads  to  a  detached  bow  shock  in  front  of  the  cone.  We  note 
that  the  solid  lines  for  large  5  3how  the  largest  variation  for  the 
different  Y’s,  while  near  the  Mach  wave  (3in  $  -  1/M^),  the  difference  is 
small.  It  is  also  true  that  the  small  angle  cones  lead  to  very  weak  oblique 
shocks  with  small  density  jumps. 


Another  way  of  presenting  the  results  on  the  maximum  3hock  angle  is  given 
in  Figure  4.  Here  the  maximum  half  angle  5  for  an  attached  shock  on  a  cone 
is  plotted  as  a  function  of  the  free  stream  Mach  number  M1  (solid  lines).  The 
corresponding  angles  for  a  two-dimensional  wedge  are  shown  as  dashed  lines. 
The  latter  results  are  simply  obtained  from  eqns.  (8)  and  (10) 
since  5  -  x  for  a  wedge.  The  large  difference  between  a  wedge  and  cone 
in  Figure  4  clearly  indicates  a  strong  dependence  on  the  geometry  of  the 
obstacle. 


Standoff  distance  for  a  detached  bow  3hock 


The  stage  is  now  set  for  computing  the  standoff  distance  for  a  detached 
bow  shock  in  front  of  a  sphere  or  a  cone  with  5  >  <5  max.  The  approximate 
procedure  is  given  by  Shapiro10  starting  on  page  884.  Since  there  are  many 
substeps  and  Shapiro's  description  is  fairly  complete  we  will  not  repeat  it 
here.  We  do  mention  several  items.  First,  the  method  requires  the 
determination  of  the  sonic  point  on  the  surface  of  the  sphere  which  depends 


upon  a  knowledge  of  5  for  a  cone  derived  in  the  last  subsection.  Second, 

max 


the  method  requires  the  shock  angle  $  for  which  the  downstream  flow  i3 

3 


exactly  sonic.  This  can  be  determined  from  eqn.  (9)  by  setting  M2  -  1  and 


solving  for  sin  $  ; 

3 


sin  <f>3  -  [ 


2YM1^(1  +  a) 


{a2  -  1  *  M^2(y  *  a2)  + 


♦  [af  -  1  ♦  M^y  ♦  a2)}2  +  4y(  1  +  a)*]1'*}]1^, 


2i1/2n1/2 


where  a  -  (Y  -  1)/2.  Third,  eqn.  (22.8)  of  Shapiro  is  not  apparent  but  can  be 
derived  using  his  notation  as  follows: 


p  V  (p  /T  )  M  C  p  T 
3  S  *S  3  S  S  _  _3  j  "  1_ 


pjm  (P.TJ  P„  Ts  M 

2 


3  ® 

Y  +  1 


.  Ik  !_  r r_S _ )  r,  .  UlJ.  -  <> 

p  M  UY  +  V  u  2 


0® 


where  we  have  used  M_  -  1 ,  the  isentropic  relations  for  a  perfect  gas, 


^  > 


Y  -  1  w2i  Y  -  1 
Pg  *  P  [1  +  —2—  M  J 


T  -  T  (l 
o 


+  ^M2), 


and  the  fact  that  the  stagnation  (subscript  o)  temperature  does  not  change 

across  a  shock.  (The  subscript  ®  denotes  free  stream  conditions,  subscript  s 

denotes  sonic  line  conditions,  and  subscript  c  denotes  centroid  streamlines 

values).  The  ratio  p  /p  can  be  obtained  by  using  eqn.  (4)  for  the  centroid 

oc  o® 

streamline  at  the  shock  front  and  then  applying  the  isentropic  relations  for  a 
perfect  gas.  Fourth,  as  Shapiro  notes,  the  method  indicates  that  the  shape  of 
the  obstacle's  nose  upstream  of  the  sonic  point  is  not  important  in 

t 

determining  the  shock  location. 

The  outcome  of  the  procedure  gives  the  ratio  of  the  standoff  distance  to 
the  height  of  the  sonic  point  on  the  body.  The  result  can  easily  be  converted 
to  more  practical  distances  through  geometrical  conditions  depending  on  the 
body  shape.  The  solution  for  a  sphere  and  a  conical  bullet  are  presented  in 
Figures  5a  and  5b,  respectively.  We  emphasize  that  the  results  for  a  cone  are 
only  applicable  if  the  cone  half  angle  6  is  larger  than  the  5  from  Figure 

OlclX 

4.  In  going  from  an  experimentally  observed  L/d  to  a  Mach  number,  the 
resulting  M1  is  quite  sensltve  to  Y  for  small  L/d.  In  going  the  other 
direction,  i.e.,  starting  from  a  given  M1  ,  the  predicted  L/d  is  not  very 
sensitive  to  Y.  It  is  also  clear  that  the  ratio  L/d  is  not  discriminating  for 
large  M1 . 

Finally,  we  note  that  comparison  of  the  results  for  the  case  Y  -  1  .4  with 
experiments  is  summarized  by  Shapiro  and  shows  good  agreement.  In  addition,  a 
natural  question  which  arises  is  how  soon  after  the  main  shock  collides  with 
the  obstacle  is  the  above  method  applicable.  We  have  employed  the  above 
method  to  the  last  two  photographs  in  Bryson  and  Gross^  which  show  the  shock- 
sphere  interaction  in  stage  three  of  Figure  1 .  We  find  an  upstream  Mach 
number  relative  to  the  sphere  of  1.5-1 .6,  while  the  actual  Mach  number  is 
slightly  less,  i.e.,  1.333.  (The  actual  value  was  obtained  from  the  quoted 
incident  Mach  number  in  the  Bryson  and  Gross  experiment  and  the  shock  jump 
relations) . 


III.  Applicability  to  the  Experiment 


As  described  in  the  Introduction,  we  are  proposing  that  a  small  obstacle 
in  the  upgraded  experiment  be  used  to  probe  the  nature  of  the  cavity.  The 
flow  Mach  number  may  then  be  estimated  from  the  standoff  distance  of  the  bow 
shock  and  Figures  5a  and  5b.  It  is  important  to  emphasize  that  the  relevent 
interaction  is  that  of  the  obstacle  with  the  flow  in  the  cavity  and  not  the 
interaction  of  the  main  shock  with  the  obstacle.  In  the  latter  case  the  Mach 
number  of  the  main  shock  i3  large  while  in  the  former  the  Mach  number  is 
reduced  since  the  cavity  temperature  is  probably  hotter  than  that  of  the 
ambient  material  and  the  gas  velocity  is  smaller  than  the  3hock  front.  In 
general,  the  smaller  the  Mach  number  the  larger  the  variation  in  the 
geometrical  configuration  of  the  shock  interaction,  and  hence  the  easier  it  is 
to  predict  and  obtain  information  on  the  flow.  We  further  suggest  that  the 
standoff  distance  of  a  bow  shock  about  a  fat  cone  or  sphere  be  studied  rather 
than  the  attached  shock  on  a  cone.  In  the  latter  case.  Figures  3a  and  3b  3how 
that  it  is  only  for  narrow  cones  that  a  large  range  in  M1  is  allowed,  but  for 
these  type  of  obstacles  the  density  jump  across  the  oblique  shock  is  weak  and 
it  may  be  difficult  to  detect  with  shadowgraphy  techniques. 

There  are,  however,  several  caveats  that  need  to  be  considered  before  the 
theory  can  be  employed.  According  to  section  II,  the  essential  assumptions 
behind  the  theoretical  development  were  that  the  upstream  flow  was  uniform  and 
that  the  flow  on  each  side  of  the  shock  was  steady  and  adiabatic. 

Clearly,  for  an  obstacle  in  the  cavity  of  the  experiment,  the  upstream 
flow  is  not  uniform  but  diverging  since  the  blast  wave  expansion  is  more 
nearly  spherical  than  planar.  Likewise  the  flow  is  not  steady.  Both  of  these 
objections  may  be  somewhat  alleviated  by  using  small  spheres  in  the  sense  that 
d/D  is  small,  where  d  is  the  obstacle's  diameter  and  D  is  the  distance  of  the 
obstacle  from  the  target.  Spheres  are  better  than  cones  for  the  uniform 
problem  since  one  need  not  worry  about  alligning  the  cone  center  line  with  the 
flow  direction.  As  long  as  the  gradient  length  scales  of  the  flow  variables 
are  large  compared  to  d  and  the  temporal  scale  is  long  compared  to  the  flow 
transit  time,  the  problem  of  unsteady  flow  is  also  manageable.  We  will  have 
more  to  say  about  rapidly  changing  flows  at  the  end  of  this  section. 
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The  adiabaticity  of  the  cavity  flow  is  also  questionable.  If  the  cavity 
temperature  Tc  is  hot,  electron  thermal  conduction  may  be  important  and  a  two 
fluid  description  could  be  more  appropriate  for  the  cavity.  This  leads 
directly  to  the  problem  that  if  Te  is  large  (or  the  cavity  density  is  low 
enough)  the  mean  free  path  in  the  cavity  could  become  larger  than  a  small 
obstacle.  An  example  of  this  case  is  that  of  the  Moon  in  the  supersonic  solar 
wind  which  shows  no  bow  shock.  It  should  be  noted  though,  that  for  the 
theoretical  and  numerical  models  presented  on  the  experiment,  either  the 
cavity  is  so  hot  that  the  flow  is  mostly  subsonic,  in  which  case  there  is  no 
bow  shock,  or  the  cavity  is  cold  and  the  adiabatic  approximation  is 
reasonable. 

There  is  a  further  subtle  problem  mentioned  by  B.  Ripin12  which  would  be 
associated  with  any  object  in  the  experimental  chamber.  When  the  laser  is 
turned  on,  the  target  is  heated  to  such  a  degree  that  very  energetic  photons 
are  emitted.  These  photons  may  ionize  and  heat  the  surface  of  an  obstacle  and 
may  lead  to  a  flow  off  the  obstacle.  This  in  turn  would  form  a  boundary  layer 
around  the  obstacle  and  change  the  nature  of  the  shock  interaction.  At  the 
present  time  there  appears  to  be  no  method  for  testing  this  or  assessing  its 
impact. 

Given  all  the  above  mentioned  problems,  and  the  requirement  of  an  assumed 
velocity  profile  to  obtain  a  temperature  from  the  observed  Mach  number,  we 
suggest  that  it  would  still  be  of  interest  to  place  an  obstacle  in  the  path  of 
the  expanding  cavity  just  to  watch  the  temporal  change  in  the  standoff 
distance  of  the  bow  shock.  As  the  cavity  expands  the  fixed  obstacle  would 
sample  deeper  and  deeper  into  the  cavity  with  time.  The  gross  evolution  could 
then  be  compared  with  the  predictions  from  theoretical  and  numerical  models. 
For  example,  we  consider  three  distinct  models  for  the  blast  wave.  First  in 
Figure  6  we  plot  the  local  Mach  number  for  the  adiabatic  Taylor-Sedov  blast 
wave  and  the  isothermal  blast  wave  as  a  function  of  r/Rg  where  Rg  is  the 
radius  of  the  main  shock.  Since  these  theoretical  solutions  are  self-similar 
the  Mach  number  has  a  self-similar  profile.  Note  that  one  would  expect  a  bow 
shock  around  an  obstacle  only  in  the  outer  10*  of  the  cavity.  Once  the 
obstacle  gets  inside  this  region  the  high  temperature  and  low  velocity  lead  to 
subspnic  flow.  For  the  second  model,  Figure  7  shows  the  local  adabatic  Mach 
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number  at  two  times  from  Stellingwerf ’s°  model  "LSNQ."  The  parameters  include 


a  25  mn  kl  target  with  a  5  Torr  background  N£  gas  and  a  laser  energy  of  23 
Joules.  The  ratio  of  specific  heat  Y  was  evaluated  from  the  specific  internal 
energy  e  and  the  relation  Y  »  1  +  p /pc.  Note  that,  contrary  to  the  self- 
similar  models,  the  Mach  number  is  large  throughout  the  cavity.  Finally,  the 
third  model  is  a  recent  numerical  simulation^  of  the  experiment  which 
includes  the  effects  of  a  non-spherical  outflow  from  the  target  and  non¬ 
equilibrium  chemistry.  Details  will  be  presented  in  a  forth  coming  NRL  memo 
report.  Figure  8  shows  the  Mach  number  at  four  different  times  for  the 
initial  conditions  of  a  4.6  ym  thick  Al  target,  situated  at  x  »  0.0  cm  ,  2.5 
Torr  background,  and  23  Joules  of  laser  energy.  The  Figure  shows  that  at 
later  times  this  third  model  resembles  the  Taylor-Sedov  blast  wave  in  that  the 
front  edge  moving  along  the  positive  x-direction  i,e.,  back  toward  the  laser, 
is  locally  supersonic.  (The  horizontal  dashed  line  marks  a  Mach  number  of 
unity.)  The  backside  has  a  quite  irregular  structure  though.  If  an  experiment 
with  a  stationary  obstacle  can  be  constructed,  it  should  be  a  simple  matter  to 
verify  which  of  the  competing  models  are  valid. 

IV .  Summary 

We  have  partially  reviewed  the  method  for  computing  the  attached  and 
detached  3hock  geometry  for  uniform,  steady,  adiabatic,  supersonic  flow  around 
a  blunt  obstacle.  In  Figures  3a  and  3b  the  angle  of  a  conical,  attached  shock 
about  a  cone  of  half  angle  5  is  plotted  against  the  upstream  Mach  number  M1 . 
The  maximum  half-angle  5  for  which  a  cone  can  support  an  attached  shock  is 

[QaX 

3hown  in  Figure  4.  if  the  conditions  for  an  attached  shock  are  not  met. 
Figures  5a  and  5b  show  the  relation  between  the  standoff  distance  of  the 
detached  bow  shock,  and  the  upstream  Mach  number  for  spheres  and  conical 
bullets,  respectively. 

In  this  memo  report  we  suggest  that  a  small  obstacle,  preferably  a 
sphere,  be  placed  in  the  path  of  the  expanding  cavity  formed  in  the  NRL 
Laser /HANE  experiment.  Then  a  measurement  of  the  standoff  distance  of  the 
resulting  bow  shock  can  be  used  with  Figure  5  to  estimate  the  local  Mach 
number  in  the  cavity.  Assuming  a  velocity  profile  for  the  cavity  material,  a 
linear  dependence  on  radius  is  typical  of  the  models,  the  proposed  diagnostic 
can  provide  a  rough  estimate  for  the  cavity  temperature.  This  physical 
quantity  is  unknown  at  the  present  time  but  is  clearly  relevent  for 


understanding  the  relationship  between  the  experiment  and  a  real  HANE. 


Given  that  the  assumptions  listed  above  in  the  theoretical  treatment  of 
the  bow  shock  geometry  may  not  be  rigorously  met  in  the  experiment,  it  would 
still  be  of  interest  to  study  the  temporal  evolution  of  the  bow  shock  standoff 
distance  and  compare  the  observations  with  predictions  from  theoretical 
models.  For  instance,  Figure  6  shows  that  the  flow  within  the  adiabatic  and 
isothermal  self-similar  blast  wave  solutions  is  only  supersonic  in  the 
outer  -  10%  of  the  cavity.  In  the  numerical  model  presented  by  Stellingwerf 
(cf.  Figure  7)  the  cavity  is  supersonic  over  most  of  its  interior  because  it 
is  cold,  while  in  the  simulation  of  Giuliani  and  Mulbrandon  (cf.  Figure  8)  the 
situation  is  close  to  that  of  the  self-similar  solutions,  at  least  on  the 
front  side.  Thus,  even  given  the  complications  listed  in  section  III  on 
applying  the  theory  to  the  experiment,  in  simply  detecting  a  bow  shock  about 
an  obstacle  one  could  distinguish  between  completely  different  models. 
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Schematic  diagram  of  the  interactin  of  a  strong  shock  wave  with 
stationary,  spherical  obstacle.  Time  increases  downward. 


The  relation  between  the  standoff  distance  of  a  bow  shock, 
characteristic  size  of  the  blunt  obstacle,  and  upstream  Mach  number; 
a)  sphere,  b)conical  bullet. 
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The  local  Mach  number  in  the  cavity  for  the  adabatic 
(T  *  1 .2  and  y  •  5/3)*  self “similar  Taylor-Sedov 
blast  wave  solution.  The  self-similar  isothermal  solution  is 
also  shown.  a3  is  the  main  shock  radius  which  follows  eqn.  (1) 
of  the  text  and  r  is  a  spherical  coordinate. 
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The  local  Mach  number  In  the  cavity  for  Stellingwerf 's  model  "L2nQ 
shown  at  two  different  times.  Here  R  is  a  spherical  radius 
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The  local  Mach  number  in  the  cavity  for  a  numerical  model  by  Giuliani 
and  Mulbrandon.  The  A!  target  initially  sits  at  x  -  0  and  expands 
both  toward  the  laser  (at  x  -  +  «)  and  backwards.  The  horizontal 
dashed  line  marks  the  Mach  number  unity  and  between  the  vertical 
lines  the  gas  is  all  Aluminum,  while  outside  of  them  the  gas  is 
Nitrogen. 


References 

.  B.H.  Ripin,  J.  Grun,  S.  Kacenjar,  E.A.  McLean,  and  J.A.  Stamper,  NRL  Memo 
Report  #5258,  "Introduction  to  the  Laser-Hane  Experiment  and  Summary  of 
Low-Pressure  Interaction  Results,"  (1984).  (AD-A138  945) 

2.  J.A.  Stamper,  B.H.  Ripin,  E.A.  McLean,  and  S.P.  Obenschain,  NRL  Memo 
Report  #5278,  "Optical  Imaging  of  a  Coupling  Region  Between  Inter- 
Streaming  Plasmas,"  (1984). 

!  B.H.  Ripin,  J.A.  Stamper,  and  E.A.  McLean,  NRL  Memo  Report  #5279.  "Blast 
Wave  Analysis  of  High-Pressure  Coupling  Shells."  (1984). 

4.  J.A.  Stamper,  "Early-Time  High-Altitude  Working  Group  Meeting,"  held  at 
NRL,  November,  1984. 

5.  L.I.  Sedov,  Similarity  and  Dimensional  Methods  in  Mechanics ,  (Academic 
Press:  New  York,  1959). 

6.  A.  Solinger,  S.  Rappaport,  and  J.  Buff,  Ap.  £,  201,  381,  "Isothermal 
Blast  Wave  Model  of  Supernova  Remnants,"  (1975). 

'.  S.F.  Gull,  Monthly  Notices  Royal  Astr.  Soc.,  161,  47,  "  A  Numerical  Model 
of  the  Structure  and  Evolution  of  Young  Supernova  Remnants,"  (1973). 

3.  R.  Stellingwerf ,  "Early-Time  High-Altitude  Working  Group  Meeting,"  held 
at  NRL,  November,  1984. 

).  A.E.  Bryson  and  R.W.F.  Gross,  J.  of  Fluid  Mechanics,  10,  1,  "Diffraction 
of  Strong  Shocks  by  Cones,  Cylinders,  and  Spheres,"  (1961). 


1.  A.H.  Shapiro,  The  Dynamics  and  Thermodynamics  of  Compressible  Fluid  Flow, 


Vols,  I  and  II,  (Ronald  Press  Co.:  New  York,  1954) 


G.I.  Taylor  and  J.W.  Macooll,  Proc.  Royal  Soc.  London  (A_) ,  1 39 ,  278,  "The 
Air  Pressure  on  a  Cone  Moving  at  High  Speeds.  -  I."  (1933) 

B.  Ripin,  private  communication. 

J.L.  Giuliani,  Jr.  and  M.  Mulbrandon,  "Early-Time  High  Altitude  Working 
Group  Meeting,"  held  at  NRL,  May,  1985 


I 


DE?.;.?7H  -  NT  OF  DEFENSE 

ASSISTANT  SEC3E7A RY  Or  DEFENSE 
COMM,  CMD ,  COM7  7  INTELL 
WASHINGTON,  D.C.  20301 

DIRECTOR 

COMMAND  CONTROL  TECHNICAL  CENTS?. 
PENTAGON  RM  BE  535 
WASHINGTON,  D.C.  20301 
Old  ATTN  C-550 
Old  ATTN  C-312  3.  MASON 

DIRECTOR 

DEFENSE  ADVANCED  RSCH  PROJ  AGENCY 
ARCHITECT  3UILDING 
1-tOO  WILSON  3L VD  . 

ARLINGTON,  VA.  22209 
0 1 CY  ATTN  NUCLEAR 

MONITORING  RESEARCH 
0 1 C Y  ATTN  STRATEGIC  TECH  OFFICE 

DEFENSE  COMMUNICATION  ENGINEER  CENTER 
1350  WIEHLE  AVENUE 
RSSTON,  VA.  22090 
0  1  C Y  ATTN  CODE  RU10 
0 1 CY  ATTN  CODE  R812 

DEFENSE  TECHNICAL  INFORMATION  CENTER 
CAMERON  STATION 
ALEXANDRIA,  VA.  2231* 

02CY 

DIRECTOR 

DEFENSE  NUCLEAR  AGENCY 
WASHINGTON,  D.C.  20305 
0 1 CY  ATTN  STVL 
0«CY  ATTN  TITL 
0 1 CY  ATTN  DOST 
03CY  ATTN  RAAE 

COMMANDER 
FIELD  COMMAND 
DEFENSE  NUCLEAR  AGENCY 
KIRTLAMQ,  AF3 ,  NM  37115 
0 1 C Y  ATTN  FC PR 


DEFENSE  NUCLEAR  AGENCY 
SAC/3NA 

3UILDINC  2C575 
KIRTLAND  AF3 ,  NM  37115 
0 1 C Y  D.C.  THORNBURG 

DIRECTOR 

INTERSERVICE  NUCLEAR  WEAPONS  SCHOOL 
KIRTLAND  AF3 ,  NM  37115 

0 1 C Y  ATTN  DOCUMENT  CONTROL 

JOINT  CHIEFS  OF  STAFF 
WASHINGTON,  D.C.  20301 

01  CY  ATTN  .J-3  WWMCCS  EVALUATION 
OFFICE 

DIRECTOR 

JOINT  STRAT  TGT  PLANNING  STAFF 
OFFUTT  AFB 
OMAHA,  MB  53113 

01 CY  ATTN  JSTPS/JLXS 
01 CY  ATTN  JPST  G.  GOETZ 

CHIEF 

LIVERMORE  DIVISION  FLO  COMMAND  DMA 
DEPARTMENT  OF  DEFENSE 
LAWRENCE  LIVERMORE  LABORATORY 
P.O.  BOX  80S 
LIVERMORE,  CA  94550 
0 1 CY  ATTN  FCPRL 

COMMANDANT 

NATO  SCHOOL  (SHAPE) 

APO  NEW  YORK  09172 

01 CY  ATTN  U.S.  DOCUMENTS  OFFICER 

UNDER  SECY  OF  DEF  FOR  RSCH  &  ENGRG 
DEPARTMENT  OF  DEFENSE 
WASHINGTON,  D.C.  20301 

0 1  CY  ATTN  STRATEGIC  Sr  SPACE 
SYSTEMS  (OS) 

COMMAND  ER/ 0 IRECTOR 
ATMOSPHERIC  SCIENCES  LABORATORY 
U.S.  ARMY  ELECTRONICS  COMMAND 
WHITE  SANDS  MISSILE  RANGE,  NM  33002 
0 1 CY  ATTN  DELAS-SO,  F.  NILES 


DIRECTOR 

BMD  ADVANCED  TECH  CTR 
HUNTSVILLE  OFFICE 
P.O.  30X  1500 
HUNTSVILLE,  AL  35807 

0 1  C  Y  ATTN  ATC-T  MELVIN  T.  CAPPS 
0 1  C  Y  ATTN  ATC-0  W.  DAVIES 
0  1  G  Y  ATTN  ATC-R  DON  RUSS 

PROGRAM  MANAGER 
BMD  PROGRAM  OFFICE 
5001  EISENHOWER  AVENUE 
ALEXANDRIA,  VA  22333 

01 CY  ATTN  DACS-BMT  J.  SHEA 

CHIEF  C-E-  SERVICES  DIVISION 
U.S.  ARMY  COMMUNICATIONS  CMD 
PENTAGON  RM  1 B269 
WASHINGTON,  D.C.  20310 

0 1 C Y  ATTN  C-  E-SERVICES  DIVISION 

COMMANDER 

FRADCOM  TECHNICAL  SUPPORT  ACTIVITY 
DEPARTMENT  OF  THE  ARMY 
FORT  MONMOUTH,  N.J.  07703 

0 1 C Y  ATTN  DRSEL-NL-RD  H.  BENNET 
0 1 C Y  ATTN  DRSEL-PL-ENV  H.  BOMKE 
0 1 C Y  ATTN  J.E.  QUIGLEY 

COMMANDER 

U.S.  ARMY  COMM-ELEC  ENGRG  INSTAL  AG Y 
FT.  HUACHUCA,  AZ  85613 

0 1 C Y  ATTN  CCC-EMEO  GEORGE  LANE 

COMMANDER 

U.S.  ARMY  FOREIGN  SCIENCE  &  TECH  CTR 
220  7TH  STREET,  NE 
CHARLOTTESVILLE,  VA  22901 
0 1 C Y  ATTN  DRXST-SD 

COMMANDER 

U.S.  ARMY  MATERIAL  DEV  &  READINESS  CMD 
5001  EISENHOWER  AVENUE 
ALEXANDRIA,  VA  22333 

0 1 C Y  ATTN  DRCLDC  J.A.  BENDER 

COMMANDER 

U.S.  ARMY  NUCLEAR  AND  CHEMICAL  AGENCY 
7500  3ACKLICK  ROAD 
BLDG  2073 

SPRINGFIELD,  VA  22150 
0 1 C  Y  ATTN  LIBRARY 


DIRECTOR 

U.S.  ARMY  BALLISTIC  RESEARCH 
LABORATORY 

ABERDEEN  PROVING  GROUND,  MD  21005 
0 1 CY  ATTN  TECH  LIBRARY, 

EDWARD  3AICY 

COMMANDER 

U.S.  ARMY  SATCOM  AGENCY 
FT.  MONMOUTH,  NJ  07703 

0 1 C Y  ATTN  DOCUMENT  CONTROL 

COMMANDER 

U.S.  ARMY  MISSILE  INTELLIGENCE  AGENCY 
REDSTONE  ARSENAL,  AL  35809 
0 1 CY  ATTN  JIM  GAMBLE 

DIRECTOR 

U.S.  ARMY  TRADOC  SYSTEMS  ANALYSIS 
ACTIVITY 

WHITE  SANDS  MISSILE  RANGE.  NM  88002 
01 CY  ATTN  ATAA-SA 
0 1 C Y  ATTN  TCC/F.  PAYAN  JR. 

0 1 C Y  ATTN  ATTA-TAC  LTC  J.  HESSE 

COMMANDER 

NAVAL  ELECTRONIC  SYSTEMS  COMMAND 
WASHINGTON,  D.C.  20360 

0 1 C Y  ATTN  NAVALEX  034  T.  HUGHES 
0 1 CY  ATTN  PME  117 
01 CY  ATTN  PME  1 1 7~T 
0 1 CY  ATTN  CODE  5011 

COMMANDING  OFFICER 
NAOL  INTELLIGENCE  SUPPORT  CTR 
4301  SUITLAND  ROAD,  BLDG.  5 
WASHINGTON,  D.C.  20390 

01  CY  ATTN  MR.  DUBBIN  STIC  12 
01 CY  ATTN  N ISC- 50 
0 1 C Y  ATTN  CODE  5404  J.  GALET 

COMMANDER 

NAVAL  OCCEAN  SYSTEMS  CENTER 
SAN  DIEGO,  CA  92152 
01CY  ATTN  J.  FERGUSON 


S'  AVAL  RESEARCH  LABORATORY 
WASHINGTON  ,  D.C.  20  375 

0  1  C  Y  ATTN  CODE  4700  S.L.  Ossakou, 
26  CYS  IF  UNCLASS 
( 0 1 C Y  IF  CLASS) 

ATTN  CODE  4780  J.D.  HUBA,  50 
CYS  IF  UNCLASS,  01CY  IF  CLASS 
01CY  ATTN  CODE  4701  I.  VITKOVITSKY 

0 1 C Y  ATTN  CODE  7500 

0 1 C Y  ATTN  CODE  7550 

0 1 C Y  ATTN  CODE  7580 

0 1 C Y  ATTN  CODE  7551 

01CY  ATTN  CODE  7555 

0 1 C Y  ATTN  CODE  4730  E.  MCLEAN 

0 1 C Y  ATTN  CODE  4108 

0 1 C Y  ATTN  CODE  4730  B.  RIPIN 

20CY  ATTN  CODE  2628 

COMMANDER 

NAVAL  SPACE  SURVEILLANCE  SYSTEM 
DAHLGREN,  VA  22448 

0 1 C Y  ATTN  CAPT  J.H.  BURTON 

OFFICER- IN-CHARGE 
NAVAL  SURFACE  WEAPONS  CENTER 
WHITE  OAK,  SILVER  SPRING,  MD  20910 
0 1 C Y  ATTN  CODE  F 3 1 

DIRECTOR 

STRATEGIC  SYSTEMS  PROJECT  OFFICE 
DEPARTMENT  OF  THE  NAVY 
WASHINGTON,  D.C.  20376 
0 1 C Y  ATTN  NSP-2141 

0 1 C Y  ATTN  NSSP-2722  FRED  WIMBERLY 
COMMANDER 

NAVAL  SURFACE  WEAPONS  CENTER 
DAHLGREN  LABORATORY 
DAHLGREN,  VA  22448 

0 1 C Y  ATTN  CODE  DF-14  R.  BUTLER 

OFFICER  OF  NAVAL  RESEARCH 
ARLINGTON,  VA  22217 
0 1 C Y  ATTN  CODE  465 
0 1 C Y  ATTN  CODE  461 
0 1 C Y  ATTN  CODE  402 
0 1 C Y  ATTN  CODE  420 
0 1 C Y  ATTN  CODE  421 


COMMANDER 

AEROSPACE  DEFENSE  COMMAND/DC 
DEPARTMENT  OF  THE  AIR  FORCE 
ENT  AF 3 ,  CO  80912 

01 CY  ATTN  DC  MR.  LONG 

COMMANDER 

AEROSPACE  DEFENSE  COMMAND/XPD 
DEPARTMENT  OF  THE  AIR  FORCE 
ENT  AF 3 ,  CO  80912 
0 1 C Y  ATTN  XPDQQ 
0 1 C Y  ATTN  XP 


AIR  FORCE  GEOPHYSICS  LABORATORY 


HANSCOM 

AFB  , 

MA 

01  731 

01  CY 

ATTN 

OPR 

HAROLD  GARDNER 

01  CY 

ATTN 

LKB 

KENNETH 

S.W.  CHAMPION 

01  CY 

ATTN 

OPR 

ALVA  T.  STAIR 

01  CY 

ATTN 

‘PHD 

JURGEN  BUCHAU 

01  CY 

ATTN 

PHD 

JOHN  P.  MULLEN 

AF  WEAPONS  LABORATORY 

KIRTLAND 

AFT 

,  NM 

8711  7 

01  CY 

ATTN 

SUL 

01  CY 

ATTN 

CA 

ARTHUR  H.  GUENTHER 

01CY 

ATTN 

NTYCE  1 LT .  G.  KRAJEI 

AFTAC 

PATRICK 

AFB, 

FL 

32925 

01  CY 

ATTN 

TN 

AIR  FORCE  AVIONICS  LABORATORY 
WRIGHT-PATTERSON  AFB ,  OH  45433 
0 1 C Y  ATTN  A AD  WADE  HUNT 
01CY  ATTN  A AD  ALLEN  JOHNSON 

DEPUTY  CHIEF  OF  STAFF 
RESEARCH,  DEVELOPMENT,  &  ACQ 
DEPARTMENT  OF  THE  AIR  FORCE 
WASHINGTON,  D.C.  20330 
0 1 C Y  ATTN  AFRDQ 

HEADQUARTERS 

ELECTRONIC  SYSTEMS  DIVISION 
DEPARTMENT  OF  THE  AIR  FORCE 
HANSCOM  AFB,  MA  01731"5000 
01 CY  ATTN  J.  DEAS 
ESD/SCD-4 


COMMAND  ER 

FOREIGN  TECHNOLOGY  DIVISION,  AFSC 
WRIGHT-PATTERSON  AFB ,  OH  *45433 
0 1 C Y  ATTN  NICD  LIBRARY 

C1CY  ATTN  ETDP  B.  BALLAP.D 

COMMANDER 

ROME  AIR  DEVELOPMENT  CENTER,  AFSC 
GRIFFISS  AFB,  NY  13441 

0 1 C Y  ATTN  DOC  LIBRARY/TSLD 
0 1 CY  ATTN  OCSE  V.  COYNE 

STRATEGIC  AIR  COMMAND/XPFS 
OFFUTT  AFB,  NB  68113 

0 1 C Y  ATTN  ADWATE  MAJ  BRUCE  BAUER 
0 1 CY  ATTN  NRT 

0 1 C Y  ATTN  DOK  CHIEF  SCIENTIST 

SAMSO/SK 
P.O.  BOX  92960 
WORLDWAY  POSTAL  CENTER 
LOS  ANGELES,  CA  90009 

0 1 C Y  ATTN  SKA  (SPACE  COMM  SYSTEMS) 
M.  CLAVIN 

SAMSO/MN 

NORTON  AFB,  CA  92409 
(M I NUT EM AN  ) 

0 1 C Y  ATTN  MNNL 

COMMANDER 

ROME  AIR  DEVELOPMENT  CENTER,  AFSC 
HANSCOM  AFB,  MA  01731 

0 1 C Y  ATTN  EEP  A.  LORENTZEN 

DEPARTMENT  OF  ENERGY 
LI3RARY  ROOM  G-042 
WASHINGTON,  D.C.  20545 

0 1 CY  ATTN  DOC  CON  FOR  A.  LABOWITZ 

DEPARTMENT  OF  ENERGY 
ALBUQUERQUE  OPERATIONS  OFFICE 
P.O.  BOX  5400 
ALBUQUERQUE,  NM  871 15 

0 1 C Y  ATTN  DOC  CON  FOR  D.  SHERWOOD 

EG&G ,  INC. 

LOS  ALAMOS  DIVISION 

P.O.  BOX  809 

LOS  ALAMOS,  NM  85544 

MCY  ATTN  DOC  CON  FOR  J.  BREEDLOVE 


UNIVERSITY  OF  CALIFORNIA 
LAWRENCE  LIVERMORE  LABORATORY 
P.O.  BOX  808 
LIVERMORE,  CA  94550 

01 CY  ATTN  DOC  CON  FOR  TECH  INFO 
DEPT 

0 1 C Y  ATTN  DOC  CON  FOR  L-389  R.  OTT 

0 1 CY  ATTN  DOC  CON  FOR  L~31  R.  HAGER 

LOS  ALAMOS  NATIONAL  LABORATORY 

P.O.  BOX  1663 

LOS  ALAMOS,  NM  87545 

0 1 CY  ATTN  DOC  CON  FOR  J.  WOLCOTT 

0 1 C Y  ATTN  DOC  CON  FOR  R.F.  TASCHEK 

0 1 C Y  ATTN  DOC  CON  FOR  E.  JONES 

01 CY  ATTN  DOC  CON  FOR  J.  MALIK 

0 1 CY  ATTN  DOC  CON  FOR  R.  JEFFRIES 

0 1 C Y  ATTN  DOC  CON  FOR  J.  ZINN 

0 1 CY  ATTN  DOC  CON  FOR  D.  WESTERVELT 

0 1 C Y  ATTN  D.  SAPPENFIELD 

SANDIA  LABORATORIES 
P.O.  BOX  5800 
ALBUQUERQUE,  NM  87115 

0 1 C Y  ATTN  DOC  CON  FOR  W.  BROWN 
01 CY  ATTN  DOC  CON  FOR  A. 

'  THORNBROUGH 

01CY  ATTN  DOC  CON  FOR  T.  WRIGHT 

0 1 C Y  ATTN  DOC  CON  FOR  D.  DAHLGREN 

0 1 C Y  ATTN  DOC  CON  FOR  3141 
0 1 C Y  ATTN  DOC  CON  FOR  SPACE  PROJECT 
DIV 

SANDIA  LABORATORIES 
LIVERMORE  LABORATORY 
P.O.  BOX  969 
LIVERMORE,  CA  94550 

0 1 C Y  ATTN  DOC  CON  FOR  B.  MURPHEY 

0 1 C Y  ATTN  DOC  CON  FOR  T.  COOK 

OFFICE  OF  MILITARY  APPLICATION 
DEPARTMENT  OF  ENERGY 
WASHINGTON,  D.C.  20545 

0 1 CY  ATTN  DOC  CON  DR.  YO  SONG 


OTHER  GOVERNMENT 


INSTITUTE  FOR  TELECOM  SCIENCES 

NAT  10 HAL 

TELECOMMUNICATIONS  & 

INFO 

ADMIN 

BOULDER  , 

CO  80303 

0  1  CY 

ATTN  A.  JEAN  (UNCLASS 

ONLY  ) 

01  CY 

ATTN  W.  UTLAUT 

01  CY 

ATTN  D.  CROMBIE 

01  CY 

ATTN  L.  BERRY 

NATIONAL  OCEANIC  &  ATMOSPHERIC  ADMIN 
ENVIRONMENTAL  RESEARCH  LABORATORIES 
DEPARTMENT  OF  COMMERCE 
BOULDER,  CO  30302 
01CY  ATTN  R.  GRUBB 
0 1  C  Y  ATTN  AERONOMY  LAB  G.  REID 


DEPARTMENT  OF  DEFENSE  CONTRACTORS 

AEROSPACE  CORPORATION 
P.O.  BOX  92957 
LOS  ANGELES,  CA  90009 
01CY  ATTN  I.  GARFUNKEL 
0 1 C Y  ATTN  T.  SALMI 
0 1 C Y  ATTN  V.  JOSEPHSON 
01 CY  ATTN  S.  BOWER 
0 1 CY  ATTN  D.  OLSEN 

ANALYTICAL  SYSTEMS  ENGINEERING  CORP 
5  OLD  CONCORD  ROAD 
BURLINGTON,  MA  01803 

0 1 C Y  ATTN  RADIO  SCIENCES 

AUSTIN  RESEARCH  ASSOC.,  INC. 

1901  RUTLAND  DRIVE 
AUSTIN,  TX  78758 

01 CY  ATTN  L.  SLOAN 
0 1 C Y  ATTN  R.  THOMPSON 

BERKELEY  RESEARCH  ASSOCIATES,  INC. 
P.O.  BOX  983 
BERKELEY,  CA  94701 

01 CY  ATTN  J.  WORKMAN 
01 CY  ATTN  C.  PRETTIE 
0 1 CY  ATTN  S.  BRECHT 


BOEING  COMPANY,  THE 
P.O.  BOX  3707 
SEATTLE,  WA  98124 

0 1 CY  ATTN  G.  KEISTER 
0 1 C Y  ATTN  D.  MURRAY 
0 1 C Y  ATTN  G.  HALL 
01 CY  ATTN  J.  KENNEY 

CHARLES  STARK  DRAPER  LABORATORY,  INC. 
555  TECHNOLOGY  SQUARE 
CAMBRIDGE,  MA  02139 
0 1  C Y  ATTN  D.B.  COX 
0 1  C Y  ATTN  J.P.  GILMORE 

COMSAT  LABORATORIES 
LINTHICUM  ROAD 
CLARKSBURG,  MD  20734 
01CY  ATTN  G.  HYDE 

CORNELL  UNIVERSITY 

DEPARTMENT  OF  ELECTRICAL  ENGINEERING 
ITHACA,  NY  14850 

0 1 C Y  ATTN  D.T.  FARLEY,  JR. 

ELECTROSPACE  SYSTEMS,  INC. 

BOX  1359 

RICHARDSON,  TX  75080 
0 1 CY  ATTN  H.  LOGSTON 
0 1 C Y  ATTN  SECURITY  (PAUL  PHILLIPS) 

EOS  TECHNOLOGIES,  INC. 

606  Wil3hire  Blvd. 

Santa  Monica,  Calif  90401 
0 1 C Y  ATTN  C.B.  GABBARD 
0 1 C Y  ATTN  R.  LELEVIER 

ESL,  INC. 

495  JAVA  DRIVE 
SUNNYVALE,  CA  94086 
01 CY  ATTN  J.  ROBERTS 
0 1 C Y  ATTN  JAMES  MARSHALL 

GENERAL  ELECTRIC  COMPANY 
SPACE  DIVISION 
VALLEY  FORGE  SPACE  CENTER 
GODDARD  BLVD  KING  OF  PRUSSIA 
P.O.  BOX  8555 
PHILADELPHIA,  PA  19101 
0 1 CY  ATTN  M.H.  BORTNER 
SPACE  SCI  LAB 


GENERAL  ELECTRIC  TECH  SERVICES 
CO.,  INC. 

HMES 

COURT  STREET 
SYRACUSE,  NY  13201 

01CY  ATTN  G.  MILLMAN 

GEOPHYSICAL  INSTITUTE 
UNIVERSITY  OF  ALASKA 
FAIRBANKS,  AK  99701 

(ALL  CLASS  ATTN:  SECURITY  OFFICER) 
01CY  ATTN  T.N.  DAVIS  (UNCLASS  ONLY) 
0 1 C Y  ATTN  TECHNICAL  LIBRARY 
0 1 C Y  ATTN  NEAL  BROWN  (UNCLASS  ONLY) 

GTE  SYLVANIA,  INC. 

ELECTRONICS  SYSTEMS  GRP-EASTERN  DIV 
77  A  STREET 
NEEDHAM,  MA  02194 

0 1 CY  ATTN  DICK  STEINHOF 

HSS,  INC. 

2  ALFRED  CIRCLE 
BEDFORD,  MA  01730 

0 1 CY  ATTN  DONALD  HANSEN 

ILLINOIS,  UNIVERSITY  OF 
107  COBLE  HALL 
150  DAVENPORT  HOUSE 
CHAMPAIGN,  IL  61820 

(ALL  CORRES  ATTN  DAN  MCCLELLAND) 

0 1 C Y  ATTN  K.  YEH 

INSTITUTE  FOR  DEFENSE  ANALYSES 
1801  NO.  BEAUREGARD  STREET 
ALEXANDRIA,  VA  22311 
0 1 CY  ATTN  J.M.  AEIN 
0 1 C Y  ATTN  ERNEST  BAUER 
01 CY  ATTN  HANS  WOLFARD 
0  1  C Y  ATTN  JOEL  BENGSTON 

INTL  TEL  4  TELEGRAPH  CORPORATION 
500  WASHINGTON  AVENUE 
NUTLEY,  NJ  07110 

0 1 C Y  ATTN  TECHNICAL  LIBRARY 

JAYCOR 

11011  TORREYANA  ROAD 

P.O.  BOX  85154 

SAN  DIEGO,  CA  921 38 

0 1  CY  ATTN  J.L.  SPERLING 


JOHNS  HOPKINS  UNIVERSITY 
APPLIED  PHYSICS  LABORATORY 
JOHNS  HOPKINS  ROAD 
LAUREL,  MD  20810 

0 1 C Y  ATTN  DOCUMENT  LIBRARIAN 
01 CY  ATTN  THOMAS  POTEMRA 
0 1 C Y  ATTN  JOHN  DASSOULAS 

KAMAN  SCIENCES  CORP 
P.O.  BOX  7463 

COLORADO  SPRINGS,  CO  80933 
0 1 CY  ATTN  T.  MEAGHER 

KAMAN  TEMPO-CENTER  FOR  ADVANCED 
STUDIES 

816  STATE  STREET  (P.O  DRAWER  QQ) 

SANTA  BARBARA,  CA  93102 
01 CY  ATTN  DASIAC 
0 1 CY  ATTN  WARREN  S.  KNAPP 
0 1 CY  ATTN  WILLIAM  MCNAMARA 
0 1 CY  ATTN  B.  .GAMBILL 

LINKA3IT  CORP 

irthe-i  orteELLS 

SAN  DIEGO,  CA  92121 

0 1 CY  ATTN  IRWIN  JACOBS 

LOCKHEED  MISSILES  4  SPACE  CO.,  INC 
P.O.  BOX  504 
SUNNYVALE,  CA  94088 
0 1 CY  ATTN  DEPT  60-12 
01 CY  ATTN  D.R.  CHURCHILL 

LOCKHEED  MISSILES  4  SPACE  CO.,  INC. 
3251  HANOVER  STREET 
PALO  ALTO,  CA  94304 

0 1  C Y  ATTN  MARTIN  WALT  DEPT  52-12 
0 1 C Y  ATTN  W.L.  IMHOF  DEPT  52-12 
0 1 CY  ATTN  RICHARD  G.  JOHNSON 
DEPT  52-12 

01CY  ATTN  J.B.  CLADIS  DEPT  52-12 

MARTIN  MARIETTA  CORP 
ORLANDO  DIVISION 
P.O.  BOX  5837 
ORLANDO,  FL  32805 

01CY  ATTN  R.  HEFFNER 


M.I.T.  LINCOLN  LABORATORY 
P.O.  BOX  73 
LEXINGTON,  MA  02173 

0 1 C Y  ATTN  DAVID  M.  TOWLE 
0 1 C  Y  ATTN  L.  LOUGHLIN 
0 1 C Y  ATTN  D.  CLARK 

MCDONNEL  DOUGLAS  CORPORATION 
5301  BOLSA  AVENUE 
HUNTINGTON  BEACH,  CA  92647 
0 1 C  Y  ATTN  N.  HARRIS 
0 1 C Y  ATTN  J.  MOULE 
0 1 C Y  ATTN  GEORGE  MROZ 
0 1 C Y  ATTN  W.  OLSON 
0 1 C Y  ATTN  R.W.  HALPRIH 
0 1 C Y  ATTN  TECHNICAL 

LIBRARY  SERVICES 

MISSION  RESEARCH  CORPORATION 
735  STATE  STREET 
SANTA  BAR3ARA,  CA  93101 
0 1 C Y  ATTN  P.  FISCHER 
0 1 C Y  ATTN  W.F.  CREVIER 
0 1 CY  ATTN  STEVEN  L.  GUTSCHE 
0 1 C Y  ATTN  R.  BOGUSCH 
0 1 CY  ATTN  R.  HENDRICK 
0 1 C Y  ATTN  RALPH  KILB 
0 1 CY  ATTN  DAVE  SOWLE 
0 1 C Y  ATTN  F.  FAJEN 
0 1 C Y  ATTN  M.  SCHEIBE 
0 1 C Y  ATTN  CONRAD  L.  LONGMIRE 
0 1 C Y  ATTN  B.  WHITE 
0 1 C Y  ATTN  R.  STAGAT 

MISSION  RESEARCH  CORP. 

1720  RANDOLPH  ROAD,  S.E. 
ALBUQUERQUE,  NEW  MEXICO  87106 
0  1  C Y  R.  STELLINGWERF 
0 1 CY  M.  ALME 
0  1  C Y  L.  WRIGHT 

MITRE  CORP 

WESTGATE  RESEARCH  PARK 
1820  DOLLY  MADISON  BLVD 
MCLEAN,  VA  22101 
01 CY  ATTN  W.  HALL 
0  1  C Y  ATTN  W.  FOSTER 

PACIFIC-SIERRA  RESEARCH  CORP 
12340  SANTA  MONICA  BLVD. 

LOS  ANGELES,  CA  90025 

0 1 C Y  ATTN  E.C.  FIELD.  JR. 


PENNSYLVANIA  STATE  UNIVERSITY 
IONOSPHERE  RESEARCH  LAB 
318  ELECTRICAL  ENGINEERING  EAST 
UNIVERSITY  PARK,  PA  16802 
(NO  CLASS  TO  THIS  ADDRESS) 

01 CY  ATTN  IONOSPHERIC  RESEARCH  LAB 

PHOTOMETRICS,  INC. 

4  ARROW  DRIVE 
WOBURN,  MA  01801 

0 1 C  Y  ATTN  IRVING  L.  KOFSKY 

PHYSICAL  DYNAMICS,  INC. 

P.O.  BOX  3027 
BELLEVUE,  WA  98009 

0 1 C Y  ATTN  E.J.  FREMOUW 

PHYSICAL  DYNAMICS,  INC. 

P.O.  BOX  10367 
OAKLAND,  CA  94610 
ATTN  A.  THOMSON 

RAD  ASSOCIATES 
P.O.  BOX  9695 
MARINA  DEL  REY ,  CA  90291 
0 1 CY  ATTN  FORREST  GILMORE 
0 1 CY  ATTN  WILLIAM  B.  WRIGHT,  JR. 

0 1 C Y  ATTN  WILLIAM  J.  KARZAS 
0 1 C Y  ATTN  H.  ORY 
0 1 CY  ATTN  C.  MACDONALD 

RAND  CORPORATION,  THE 
1700  MAIN  STREET 
SANTA  MONICA,  CA  90406 
01 CY  ATTN  CULLEN  CRAIN 
0 1 C Y  ATTN  ED  BEDROZI AN 

RAYTHEON  CO. 

528  BOSTON  POST  ROAD 
SUDBURY,  MA  01776 

0 1 CY  ATTN  BARBARA  ADAMS 

RIVERSIDE  RESEARCH  INSTITUTE 
330  WEST  42nd  STREET 
NEW  YORK,  NY  10036 

0 1 C Y  ATTN  VINCE  TRAPANI 


SC IENCE 

APPLICATIONS,  INC. 

TRW  DEFE 

N5E  h.  SPACE  SYS  CROUP 

1150  PROSPECT  PLAZA 

ONE  SPAC 

E  PARK 

LA  JOLLA 

,  CA 

92037 

R EDCNDO 

BEACH,  CA  9  0  2  ~  8 

01  CY 

ATTN 

LEWIS  M.  LINSON 

0  1  C  Y 

ATTN  R.  K.  PLSE’JCH 

01  CY 

ATTN 

DANIEL  A.  HAMLIN 

0  1  CY 

ATTN  S.  ALTSCHULER 

01  CY 

ATTN 

E.  FRIEMAN 

01  CY 

ATTN  D.  D  E- 

0  1  C  Y 

ATTN 

E.A.  STRAKER 

0  1  CY 

ATTN  D/  STOCKWLLL 

01  CY 

ATTN 

CURTIS  A.  SMITH 

SNTF / 1575 

SCIENCE 

APPLICATIONS,  INC 

VISIDYNE 

1710  GOODRIDGE  DR. 

SOUTH  3E 

D  FORD  STREET 

MCLEAN  , 

V  A 

221  02 

BURLINGTON,  MASS  01803 

0  1  C  Y 

J.  COCKAYNE 

01  CY 

ATTN  W.  REIDY 

01  CY 

E.  HYMAN 

01  CY 

ATTN  J.  CARPENTER 

01  CY 

ATTN  C.  HUMPHREY 

SRI  INTERNATIONAL 
333  RAVENSWOOD  AVENUE 
MENLO  PARK,  CA  94025 


ATTN  J.  CASPER 
ATTN  DONALD  NEILSON 
ATTN  ALAN  BURNS 
ATTN  G.  SMITH 
ATTN  R.  TSUNODA 
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